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1. Summary : It sometimes happens that a physical phenomenon (target) 



leaves occasional traces of its existence as patches of emitted material 
(blobs) that may temporarily grow in size, but which ultimately dissipate 
and become undetectable. In this report we investigate the chances of 
contacting one of these blobs, assuming that the target moves in its medium 
enough to prevent effective overlap of one blob on another. Search is assumed 
to be random within some confined region by a sensor capable of detecting a 
given concentration of emitted material; the effect of this assumption is to 
make the contact rate at any time proportional to the total projected length 
(in two dimensions) or area (in three dimensions) of all extant blobs. The 
random search assumption is made for analytical convenience; nonetheless, 
the assumption is often a surprisingly accurate representation of real 
searches, particularly when the target is moving. 

Let C be the mean number of times that a blob is contacted. Formulas 
for C are derived in sections 3 and 4 for the cases of diffusion in two and 
three dimensions. As long as C is much smaller than 1, the mean time E(T) 
until the first blob is contacted is approximately an exponential random 
variable with mean 1/AC, where X is the Poisson rate at which blobs are 
emitted. Actually, 1/XC is a slight underestimate of E(T), as is shown 
in sec. 2. Exact formulas for the probability distribution of T (formula 
A4) and E(T) (formula A5) valid for all C are derived in the appendix. 

The blob that is detected may be an old one if blobs have a long 

lifetime. The time required to detect a blob whose age does not exceed t^ 

could be analyzed by truncating the rate of detection function at 

t : l.e., r(x) =0 for x > t . The effect would be to reduce the 

o o 

dimensionless constant in formulas (9) and (16). 



2 . The mean time to detection exceeds 1/XC 



On account of the random search assumption, the total number of contacts 
when the blob emission times are given is a non-homo geneous Poisson process. 
Let Y(t) be the mean value function of this process, and let X(t) be 
the actual number of contacts up to time t. Since Y(t) is itself a 
stochastic process, X(t) is doubly stochastic [1]. Since X(t) is a 
Poisson random variable with mean Y(t), E(x(t)|Y(t)) == Y(t), and 

E(x(t)) = E(Y(t)) by conditional expectations. This is handy because 
E(Y(t)l is easy to predict in equilibrium. In fact, let C be the average 
number of contacts per blob (see sections 3 and 4 for examples of how C 
can be computed). Assuming that blobs are emitted at the Poisson rate X, 
we must have 



(1) E(x(t)) = E(Y(t)) = XCt 



It follows that the average time between contacts is 1/XC. However, this 
is not the same thing as the mean time E(T) until the first contact, 
since the origin of T is meant to be an arbitrary point in the process, 
rather than a point at which a detection has just occurred and for which 
Y(t) is therefore likely to be large. 

00 

Since T is a positive random variable, E(T) = / P(T > t)dt, and 
" 0 

therefore E(T) = E(/ P(T > t|Y(t))dt) by conditional expectations. But 

0 

the event (T > t) is the same as the event (X(t) = 0) , and 
P(X(t) = 0|Y(t)) = exp(-Y(t)), so 



( 2 ) 



E(T) = E 



/ exp(-Y(t))dt >_ / exp(-E(Y(t)))dt 
0 0 



= / exp(-XCt)dt = 1/AC , 
0 



2 



where the inequality is Jensen’s inequality applied to the convex function 
exp(-x). So the mean time to the first contact exceeds 1/AC, in general. 

The condition for E(T) to be approximately 1/AC is that Y(t) should 
not be very random’. Intuitively, this will be the case if the contribution 
due to a single blob is small; i.e., if C « 1. This is correct, but the 
proof is not trivial, so we defer it to the appendix. 



3 . Evaluation of C for turbulent diffusion in the ocean 

In [4], it is stated that the diffusion of substance from a point in the 
ocean is essentially a two dimensional problem, since there is very little 
vertical mixing compared to turbulent horizontal processes. Several radially 
symmetric formulas are offered to describe the way in which the concentration 
of a substance s(t,r) changes with distance (r) and time (t) from the 
point of injection, most of which are special cases of 



(3) 



s(t,r) = 



M 

TT(3t)^°‘^’^r(i+2/b) 



exp [-r*^/ , 



where M is the amount of material released, a and b are dimensionless, 

and 3 has whatever units are required to make s(t,r) have units of 

material per unit area. F is the Gamma function, and is required in (3) 

in order to make the volume underneath s(t,r) equal to M for all t. 

The case of classical diffusion is a = 1, b = 2. 

At time t after injection, a concentration sensor will be able to 

detect the substance within some distance p(t) of the point at which it was 

released. Suppose that M units of material could just barely be detected 

inside a circle of radius rQ if the material were spread uniformly; r^ 

is a measure of the sensor’s sensitivity. Then p(t) must satisfy 
2 

s(t,p(t)) = M/ (TTrQ) at all times t. Using (3) to solve for p(t), we 
obtain 



3 



1/b 



(4) 



p(t) = Zn 



[(6t)^‘^^^r(l+2/b) 



except that p(t) is taken to be 0 rather than negative. If the sensor 
is conducting a random search at speed v in an area of size A, then the 
rate of contacting the blob at time t is 2vp(t)/A, since 2p(t) is 
the projection of the circular blob normal to the searcher's track. The 
average number of detections per blob is therefore ** 



(5) C = / (2vp(t)/A)dt , 

0 

Let K = rQ/r(l+2/b) , 

substitute p(t) from (4) in (5), and then substitute u = (St)^°^^^/K. 
The result is 



(6) C = (2v/A) J [Ku]^^^[Jln 

0 



Kb du 



2aB(Ku) 



l-b/2a 



where the upper limit of tu is the largest number for which the integrand is 
positive. Factoring out constants and collecting powers of Ku, (6) is 



The definite integral has a closed form expression [ 3] . Substituting this in 

(7) , we obtain 

(8) C = b r(l+l/b)/[AaB (b/2a+l/2)^'^^''^] . 

Recalling the definition of K, we finally have 



4 



(9) 



C = 



V br(l+l/b) 



Aa6 r(l+2/b)’^''^‘^'^^^^ (b/2a+l/2)^'^^^^ 



In classical diffusion, since F(3/2) = /i:/2 and F(2) = 1, (9) re- 
duces to 



( 10 ) 



c = 



V Tq 



A6 (3/2) 



3/2 • 



However, classical diffusion is not a good description of turbulent 

spreading in the ocean. According to reference [4], better descriptions 

would be a=2, b=2, 6=2 cm. /sec. (Okubo and Pritchard), or a=3, b=2, 

2/3 

6=-05 cm /sec, (Okubo). Both of these models are from table 5 of 

4 3 10 2 

reference [4]. If v=10 cm/sec., ^^=10 cm., and A=10 cm , (9) gives 

C = .44 for the first of these models, and C = 1.55 for the second. These 

calculations may give the reader some idea of circumstances where C « 1. 
The fact that different investigators have arrived at models that provide 
answers differing by a factor of 3.5 should not be disturbing when the sub- 
ject is turbulence in the ocean. 

Figure 1 shows (4) for the same two models. The maximum blob size 
depends only on r^ and b, and is consequently the same in both cases. The 

fact that the area under one curve is 3.5 times as large as the area under 

the other is essentially a question of blob lifetime. 



4. Classical diffusion in three dimensions . 

This derivation will parallel that of sec. 3. If S(t,r) is now the 
concentration of material per unit volume and D the diffusion constant, 
then (3) is replaced by 

(11) S(t,r) = ^ exp(-r^/2Dt) . 

(2TTDt)^' 



We also have 



5 
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( 12 ) 



1 



(4/3)7Tr^ . 

:^exp(-p (t)/2Dt), 

(277Dt)-^' 

where is the radius of the largest sphere within which M units of 

material would be detectable if spread evenly, 

2 1/3 

Solving (12) for p (t), and defining K = (9 tt/ 2) , we have 



(13) 



p (t) = 2Dt Jin 



1 3/2 
0 



KDt 



= 3Dt In 



KDt] 



2 

except that p(t) = 0 for t > r^/KD, The rate of detection at time t 
is now r(t) = v(7Tp (t))/V, where V is the volume of the confining region, 
so 

00 

(14) C = TTv/V / p (t) dt. 

0 

2 

Substituting u = KOt/r^, 

(15) C = (7Tv/V)(r^/KD)(3r^/K) / u iln(l/u) du 

u ^0 

The definite integral [3] is 1/4, so 



(16) 



C = .1343 (vr^/VD). 
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Appendix : 

Let r(x) be the detection rate on a blob at a time x after it has 
been emitted, and suppose blobs are emitted in a Poisson process with rate X. 
Let N(t) be the number of blobs emitted in [0,t], and let be the times 

at which they are emitted, i = 1, ..., N(t). For 0 ;< u £ t, the total 

N(t) 

detection rate at time u is ^ r(u-S.), where r(x) is taken to be 0 

i=l ^ 

for X < 0. Therefore the average number of detections in the interval 
[t-A,t], given the emission times, is 

N(t) t 

(Al) Y(t) = ^ / r(u-S.) du, A_> 0. 

i=l t-A ^ 

The stochastic process Y(t) is a filtered Poisson process with impulse 

t t-T 

response h(t,x) = / r(u-T) du = / r(v) dv. The generating function 

t-A t-T-A 

of Y(t) is therefore known [2]: 



f t 

(A2) = E(exp(SY(t)) = exp X j [exp(Sh(t,x)) - 1] dx 

0 



By substituting x = t - X, this can be written 



(A3) 



= exp 



X / [exp(S / r(v) dv) -1] dx 
0 x-A 



The number of detections in the interval [t-A,t] is a Poisson 
random variable with mean Y(t). The probability of no detections in this 
interval is therefore E(exp (-Y (t) ) ) = ti)^(-l). If T is the time from 
t-A to the next detection, we therefore have P(T > A) = ijj^(-l). Letting 
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^ oo 



to obtain equillibrium. 



(A4) 



P(T > A) = exp 



X 



A / [exp(- / r(v) dv) -1] dx 
0 x-A 



Since T is a positive random variable, we also have 



(A5) E(T) = / P(T > A) dA, 

0 



and this provides an explicit formula for E(T) that could be evaluated by 
numerical integration in any specific case. However, suppose the mean 



number of contacts per blob (C) is <<1. 

X 

that / r(v)dv <<1, and therefore that 
x-A 



00 

Since C = / r(v) dv, it follows 

0 



1 - exp 



- / r(v) dv f r(v) dv, and therefore that 

x-A ^ x-A 



(A6) 



P(T > A) exp(-X f(A)), where 



CO X 

(A7) f(A) = / [ / r(v) dv] dx 

0 x-A 



Since f(0) = 0, and since (d/dA) f (A) = / r(x-A) dx = C, it must be true 

0 

that f(A) = CA, and therefore that T is approximately exponential with 
mean 1/XC. Furthermore since all of the above approximations can be 
written as inequalities, 

(A8) E(T) > 1/XC, 

with equality holding when C << 1. 
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